Dynamics of magnetic charges in artificial spin ice 
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Artificial spin ice has been recently implemented in two-dimensional arrays of mesoscopic magnetic 
wires. We propose a theoretical model of magnetization dynamics in artificial spin ice under the 
action of an applied magnetic field. Magnetization reversal is mediated by domain walls carrying 
two units of magnetic charge. They are emitted by lattice junctions when the the local field exceeds 
a critical value He required to pull apart magnetic charges of opposite sign. Positive feedback from 
Coulomb interactions between magnetic charges induces avalanches in magnetization reversal. 



Spin ice fl!] shares some remarkable properties with 
water ice 2]: both possess a very large number of low- 
energy, nearly degenerate configurations satisfying the 
Bernal-Fowler ice rules. In water ice, an ion has 
two protons nearby and two farther away; in spin ice, 
two spins point into and two away from the center of ev- 
ery tetrahedron of magnetic ions. Because the ice rules 
are satisfied by a large fraction of states, the system re- 
tains much entropy down to very low temperatures 
Low-frequency dynamics in ice is associated with the mo- 
tion of defects violating the ice rules. In water ice, these 
defects carry fractional electric charges of ±0.62e (ionic 
defects) and ±0.38e (Bjerrum defects) Fractionaliza- 
tion takes an even more surprising form in spin ice: while 
the original degrees of freedom are magnetic rfipoles, the 
defects are magnetic monopoles [3-[§|- 

The charge of an ice defect is defined in terms of the 
net flux of electric field E or magnetic field H emerging 
from the defect. On the atomic scale, the flux is obscured 
by the fields of background ionic charges or magnetic 
dipoles. Coarse graining is required to reveal the field 
flux of a defect on longer length scales Q . An alternative 
approach is to alter the model by stretching point-like 
spin dipoles into dumbbell magnets until they touch one 
another, while keeping their dipole moments fixed [^. At 
the expense of a slight change in the Hamiltonian, the 
magnetic charge of a defect becomes well defined even on 
the microscopic scale. It equals ±2g = ±2/j,/a, where /i 
is the dipole moment and a is the length of a dumbbell. 

The dumbbell model is realized in artificial spin ice, a 
network of submicron ferromagnetic islands Q or wires 
lol- 1^. Each element represents a spin whose mag- 
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FIG. 1: (a) A configuration of square spin ice with no mag- 
netic charges, (b) Honeycomb spin ice always has magnetic 
charges, (c) Magnetized honeycomb spin ice. 



netic dipole moment is aligned with the wire by shape 
anisotropy. Fig. [1] The magnetostatic energy is a posi- 
tive definite quantity E^i-p = (I/Stt) / H^dV, where the 
integral is taken over the entire space. It is minimized 
when the magnetic field H = 0. -Edip can be expressed as 
the Coulomb interaction of magnetic charges with den- 
sity p(r) = V-H/47r - V-M. The field is zero, and the 
energy is minimized, when there are no magnetic charges. 
This yields the ice rule: a network node with zero mag- 
netic charge has zero influx of magnetization. The zero- 
flux rule can be satisfied in square ice. Fig. (TJa), but not 
in honeycomb ice. Fig. [Ifb), also known as kagome ice, 
where the allowed values of magnetic charge Q on a site 
are ±g and ±3g in units oi q = MA, where M is the 
magnetization of the magnetic wire and A is its cross 
section. Minimization of magnetic charge restricts Q to 
the values of ±q, yielding the modified ice rule for this 
lattice: two arrows in and one out, or vice versa 13, 



The presence of residual magnetic charges in honey- 
comb ice even at low temperatures may result in a se- 
quence of two phase transitions as its temperature is 
lowered: magnetic charge order appears first, spin order 
arises later l^, lil- Unfortunately, thermal fluctuations 
are virtually absent in artificial spin ice: reversing the 
direction of magnetization in a single wire requires going 
over an energy barrier of a few million kelvins [9] . Left to 
itself, the system remains forever in the same magnetic 
microstate. Wang et al. suggested a way to introduce 
magnetization dynamics into artificial spin ice by plac- 
ing the system in a rotating magnetic field of an oscillat- 



ing magnitude 15, T^, the analog of fiuidizing granular 
matter through vibration. It has been suggested [13, [3| 
that such induced dynamics of magnetization effectively 
create a thermal ensemble with an effective temperature. 

In this Letter we present an entirely different approach 
to the dynamics of artificial spin ice that incorporates 
the physics of magnetization reversal in ferromagnetic 
nanowires, a process mediated by the creation, propaga- 
tion, and annihilation of ma gne tic domain walls IS, 2(|. 



It is inherently dissipative [21l . l22j : as a domain wall 



propagates, magnetic energy is transferred to the lattice. 
Like fluidized granular matter, artificial spin ice is a sys- 
tem far out of equilibrium and it is not obvious that it 
can be described in the framework of equilibrium ther- 
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modynamics f23| . Mesoscopic degrees of freedom of spin 
ice tend to move downhill in the energy landscape until 
they come to rest at a local energy minimum. We use 
this approach to describe the dynamics of magnetization 
observed in honeycomb spin ice 11[ in an applied field. 

In static equilibrium, artificial spin ice is fully de- 
scribed by specifying the direction of the magnetization 
vector in every link of the lattice. These are Ising vari- 
ables because magnetization is aligned with the wire. 
Sites of the lattice carry magnetic charge of ±g or ±3q as 
explained above. Site charges can be deduced from mag- 
netization variables because the magnetic charge equals 
the net influx of magnetization. The converse is not true 
because the number of links exceeds the number of sites 
by a factor of 3/2, so the magnetic state of artificial spin 



ice cannot be described in terms of charges alone 
Spin variables must be specified for a complete descrip- 
tion. 

Transitions between static states, triggered by the ap- 
plication of an external magnetic field, involve intermedi- 
ate states in which the magnetization of one or more links 
is being reversed. At the mesoscopic level of our theory, 
such links are pictured as having two sections uniformly 
magnetized in opposite directions separated by a domain 
wall of magnetic charge Q = ±2g 21]. The reversal of 
magnetization in a link begins with the creation of a do- 
main wall at one of the link ends. The process conserves 
magnetic charge: when a site with magnetic charge —q 
emits a domain wall of charge —2q, the charge of the site 
changes to -l-g, Fig. [51 The Zeeman force — 2g/xoi? then 
pushes the domain wall to the opposite end of the link. 

The critical field required to initiate the reversal can 
be estimated as follows. A site of charge +q and a do- 
main wall —2q attract each other with a Coulomb force 
F ^ /io2q^/(47rr^) at distances r exceeding the charac- 
teristic size of the charges a. The attraction weakens for 
short distances r ^ a when the two charges merge. The 
maximum attraction is thus -Fmax ~ /^o29^/(47ra^). To 
pull the charges apart, the Zeeman force 2qB from the 
applied field must exceed i^max, giving the critical field 
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Domain walls in nanowires of submicron width w have 
the characteristic size a w O.Qw 24]. For the permalloy 
honeycomb network of Qi et al. 11] with magnetization 
M = 8.6 X 10^ A/m, width w — 110 nni and thickness 
t ~ 23 nm, fJ-oHc ~ 50 mT. 

When the magnetic field is applied at an angle 6* to a 
link, the Zeeman force comes from the longitudinal com- 
ponent Hcosd. For this reason we expected the reversal 
to occur at a higher field H{9) ~ He/ cos6. A similar 
angular dependence has been observed in magnetic wires 
with submicron width [2^ . 

To test this phenomenological model, we performed nu- 
merical simulations of magnetization reversal in a single 
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FIG. 2: Magnetization reversal in honeycomb spin ice. (a- 
d) A domain wall is emitted at one end of a link, travels to 
the other end, and gets absorbed at the junction, (e-f) If 
the applied field is sufficiently strong, a new domain wall can 
be emitted into an adjacent link triggering its magnetization 
reversal, (g-k) When a domain wall encounters a site with 
like magnetic charge, it induces the emission of a new domain 
wall into an adjacent link. 
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FIG. 3; The reversal field H of two out of three magnetic 
wires forming a junction vs. the angle between the field and 
the axis of the wire whose magnetization is being reversed. 
The lines are fits to Eq. Q with fioHc = 52.0 and 55.3 mT. 



junction of three ferromagnetic nanowires using micro- 
magnetics software package OOMMF [i^ with the cell 
size of 2 nm x 2 nm x 23 nm. The dependence H{6) is 
not symmetric. Fig. |3l and is fit well by the function 



H{9) = HJ cos{9 + a), 



(2) 



where the offset a = 19° reflects an asymmetric distribu- 
tion of ma gne tization at the junction, as we will discuss 
elsewhere [27[. The critical-fleld parameter He varied 
slightly between links reflecting small random variations 
of the width caused by lattice discretization. Two links 
of the same junction exhibited slightly different critical 
fields i?c, Fig. El 

We use these phenomenological considerations and mi- 
cromagnetics simulations to build a discrete mesoscopic 
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model of magnetization dynamics in artificial spin ice. 
We start with a fully magnetized state in which links of 
the same orientation have the same direction of magne- 
tization and magnetic charges form a staggered pattern. 
Such a state can be obtained by placing the system in a 
strong magnetic field, Fig.[Hc). In this state, each mag- 
netic wire has uniform magnetization pointing along the 
wire's axis and each junction contains a magnetic charge 
of ±1 in the units of g = Mtw determined by the flux 
of magnetization into the junction. The external field is 
then applied in the opposite direction with a gradually in- 
creasing magnitude. Magnetization reversal begins when 
the net field Hnot at one of the junctions exceeds a crit- 
ical value determined by Eq. [21 The net magnetic field 
Hnot is a superposition of the applied field Happ and of 
the demagnetizing field of the sample Hdcm- The latter 
is computed as a sum of Coulombic fields of individual 
junctions, H = Qr/{4TTr'^). The junction, initially con- 
taining charge ±1, emits a domain wall with charge ±2 
and changes its own charge to =Fl. The emitted domain 
wall is pushed by the magnetic field to the other end 
of the link, reversing the link magnetization in the pro- 
cess. Fig. [2] (b-c). Quenched disorder, inevitably present 
in real samples, is modeled by setting at random slightly 
different critical fields He in individual wires with a mean 
He and a distribution width AHc- 

As the domain wall with charge ±2 reaches the other 
end of the link, its further fate depends on sign of the 
magnetic charge it meets at the junction. If the charge 
is of opposite sign, =f1, then the domain wall is absorbed 
by the junction. Fig. [^l^c-d), whose charge reverts to ±1. 
If the net field is strong enough to stimulate the emission 
of a new domain wall of charge ±2 out of this junction. 
Fig. Efe), one of the adjacent links reverses its magneti- 
zation. Fig. [2l[f). Otherwise, the evolution stops at the 
stage shown in Fig. HJd). 

Alternatively, if the domain wall comes to a junction 
with the same sign of charge. Fig. IDJa-b), it stops dis- 
tance a short of the junction thanks to magnetostatic 
repulsion. While this could be a new equilibrium posi- 
tion, the charged domain wall creates a field of strength 
2Hc at the junction, so that the net field at the junc- 
tion is close to SHc- Its projection onto an adjacent link, 
l.5Hc, is sufficient to stimulate the emission of a new 
domain wall of charge ±2 into that link. Fig. lU^c). The 
junction, now carrying charge of the opposite sign, =f1, 
pulls in the original domain wall and settles down in a 
state with charge ±1, Fig. IH^d). 

The sequence illustrated in Fig. |4] explains why ice rule 
violations are hard to find in honeycomb ice of Qi et al. 
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FIG. 4: Magnetization reversal in uniformly magnetized spin 
ice. (a-b) In the bulk, the reversal in a link magnetized against 
the field would lead to the formation of triple charge, which 
can only happen when the field is of order 3Hc- (c) Instead, 
the reversal occurs first in links magnetized at 120° to the 
field when H ~ 2Hc- (d-g) At the edge, the reversal begins 
when H ~ He and propagates into the bulk. 



that states with charge ±3 are only transients. Much 
stronger disorder exists in cobalt samples of Ladak et 
al 
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Unless variations of the critical field are so strong 
that He at some junctions exceeds l.bHc, triply charged 
junctions, Fig.|4l^b), are unstable and decay via the emis- 
sion of a new domain wall. Fig. IDJc-d) . Permalloy sam- 
ples of Qi et al. exhibit a Gaussian distribution of critical 
fields with a standard deviation AHc = O.QAHc [28], so 



12l | who observed magnetization reversal in a field 
range between H = 50 and 75 mT. Thus some of the 
domain walls encounter junctions whose critical field ex- 
ceeds 1.5H, which explains the presence of charges ±3. 

In the limit of weak disorder, AHc ^ He, there is 
another characteristic scale of the field that becomes im- 
portant. The new scale set by the demagnetizing field 
of the sample Hdem, is the strength of the field created 
by a unit magnetic charge, Q — Mtw, at a neighboring 
junction distance L away, Hq = Mtw/ {Air L'^). When 
AHc ^ Hq, the reversal of magnetization is controlled 
mostly by the effects of quenched disorder, with links re- 
versing in a largely independent fashion in the order of in- 
creasing critical field He- Conversely, when AHc ^ ^^o, 
the reversal proceeds in a correlated fashion because of 
a positive feedback: the reversal of magnetization in one 
link redistributes magnetic charges at its ends, which in 
turn increases the net field at adjacent junctions and thus 
triggers the emission of domain walls there. In samples 
of Qi et al, Hq = 0.87 mT, which is comparable to the 
width of their reversal region, AHc = 2 mT. 

We simulated magnetization reversal in this model 
with the critical fields uniformly distributed in an in- 
terval of width AHc = 5 mT around the mean He — 50 
mT and the Coulomb field scale Hq = 0.87 mT. For sim- 
plicity we set the offset angle a = 0. A sample contain- 
ing 937 links was initially magnetized along one subset 
of links. Fig. [H^c). Subsequently, the field was switched 
off and a reversal curve M(iJ) was measured in field ro- 
tated through angle 9 from the initial direction. For 120°, 
quenched disorder dominates so that magnetization re- 
versals occur largely independently, in two stages. Links 
magnetized against the field switch when the applied field 
is within the range He ± AHc/2, whereas links magne- 
tized at 120° to the field switch in the range 2Hc ± AHc- 
The net magnetization grows in an approximately 
linear fashion in both ranges. Fig. [SI as expected for 
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We have presented a discrete model of artificial spin ice 
where magnetization dynamics is mediated by domain 
walls carrying magnetic charge. Interactions between 
magnetic charges compete with the effects of quenched 
disorder. In samples with low disorder, positive feed- 
back from charge redistribution is responsible for mag- 
netic avalanches that have been observed in some exper- 
imental situations. 

We thank John Cumings for numerous discussions and 
for sharing his unpublished data. This work was sup- 
ported in part by the NSF Grant No. DMR-0520491. 



FIG. 5: Simulated magnetization reversals. A sample is ini- 
tially magnetized in a strong field directed as in Fig. [ijc). 
Subsequently, the field is switched off and reapplied at angles 
6 = 120° and 180° to the initial direction. Vertical dashed 
lines are at He ± AHc/2 and 2Hc ± AHc- Inset: Distribution 
of avalanche lengths D{s) in the range of fields near He- De- 
viations near the bottom of the graph are due to statistical 



links with a uniform distribution of He- Links do not 
reverse completely independently from one another: as 
noted previously, the redistribution of magnetic charges 
induced by the reversal of magnetization in one link may 
trigger another reversal nearby. We observed that re- 
versals often involves small groups of links. As can be 
seen in the inset of Fig. O the distribution of the num- 
ber of links s reversing in a single event is Gaussian, 
D{s) (X exp(-sV2C^), with ^ = 4.6. 

An entirely different process is observed when the field 
is rotated through 9 = 180°. In this case. Coulomb in- 
teractions play a major role and the reversal proceeds 
through avalanches evidenced by steps in Mx{H), Fig.[5l 
When the field is near the reversal cannot begin in 
the bulk because links parallel to the applied field have 
the wrong sign of magnetic charges at the ends and will 
reverse only in a much higher field (of order 3 He). Links 
at the edges have no such problem and the reversal starts 
when a site at the edge emits a domain wall. Fig. IDJd-e) . 
When the domain wall reaches the other end of the link, 
it encounters a site with like magnetic charge and trig- 
gers the emission of a new domain wall, Fig.|4l^f), and the 
reversal of magnetization in an adjacent link. Fig. lU^g). 
This triggers an avalanche of reversals that stops when 
the traveling domain wall is absorbed by a junction with 
a large critical field He or runs into already reversed links 
2^. The distribution of avalanche lengths (Fig. [5]) 
fits a power law, D(s) oc s""^, with the exp onent t — 1.6, 
indicative of self-organized criticality [30|. Chain rever- 
sals involving 3 links have been reported by Ladak et al. 
[l3 | in this geometry; avalanches involving up to 39 links 
have been observed by Daunheimer et al. [28|. 
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